Weakly Nonlinear Theory of Dynamic Fracture 
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The common approach to crack dynamics, linear elastic fracture mechanics (LEFM), assumes 
infinitesimal strains and predicts a r~^''^ strain divergence at a crack tip. We extend this framework 
by deriving a weakly nonlinear fracture mechanics theory incorporating the leading nonlinear elastic 
corrections that must occur at high strains. This yields strain contributions "more-divergent" than 
r~^'^ at a finite distance from the tip and logarithmic corrections to the parabolic crack tip open- 
ing displacement. In addition, a dynamic length-scale, associated with the nonlinear elastic zone, 
emerges naturally. The theory provides excellent agreement with recent near-tip measurements that 
can not be described in the LEFM framework. 



Understanding the dynamics of rapid cracks is a major 
challenge in condensed matter physics. For example, high 
velocity crack tip instabilities [l|, y] remain poorly under- 
stood from a fundamental point of view. Much of our 
understanding of how materials fail stems from Linear 
Elastic Fracture Mechanics (LEFM) [3|, which assumes 
that materials are linearly elastic outside of a small zone 
where all nonlinear and dissipative processes occur. A 
central facet of LEFM is that strains diverge as r~^/^ 
at a crack's tip and that this singularity dominates all 
other strain contributions in this region. Linear elastic- 
ity should be expected to break down before dissipative 
processes occur. The small size and rapid propagation 
velocity of the near-tip region of brittle cracks have, how- 
ever, rendered quantitative measurements of the near-tip 
fields elusive. 

In the companion Letter [J] such direct near-tip mea- 
surements of the displacement field u{r) were achieved 
for Mode I cracks propagating at rapid velocities, v. 
Defining (r, 9) as coordinates moving with the crack tip, 
the propagation direction, x is defined hy 9 = and 
the loading direction, y, hy 9 = tt/2. As predicted by 
LEFM, these experiments revealed that the crack tip 
opening profile, Uy{r, ±7r), is parabolic beyond a velocity- 
dependent length-scale S{v). However, it was shown that 
although Ux {r, 9 = 0) in this range also follows the func- 
tional form predicted by LEFM, its parameters are in- 
consistent with those described by Uy{r, ±7r)! Moreover, 
the strain component ej,j,(r, 0) = dyUy{r,0) was wholly 
incompatible with LEFM, indicating a "more-divergent" 
behavior than r~^/^. These puzzling discrepancies be- 
come increasingly severe as v increases. 

In this Letter, we show that all of these puzzles can 
be quantitatively resolved by taking into account nonlin- 
ear corrections to linear elasticity, which must be relevant 
near the crack tip. This is achieved by perturbatively ex- 
panding the momentum balance equation for an elastic 
medium up to second order nonlinearities in the displace- 
ment gradients. The resulting theory provides a novel 
picture of the structure of the fields surrounding a crack 
tip, and may have implications for our understanding of 
crack dynamics. 



Nonlinear material response at the large strains near a 
crack's tip motivates us to formulate a nonlinear elastic 
dynamic fracture problem under plane stress conditions. 
Consider the deformation field (p, which is assumed to 
be a continuous, differentiable and invertible mapping 
between a reference configuration x and a deformed con- 
figuration x' such that x' = 4>{x) = x + u{x). The 
deformation gradient tensor F is defined as F = V4> or 



explicitly Fij =Si 



-djUi 



The first Fiola-Kirchhoff stress 



tensor s, that is work-conjugate to the deformation gra- 
dient F, is given as s = dFU{F), where U{F) is the strain 
energy in the deformed configuration per unit volume in 
the reference configuration [5|. The momentum balance 
equation is 



V • s = pdtt4> 



(1) 



where p is the mass density. Under steady-state prop- 
agation conditions we expect all of the fields to de- 
pend on X and t through the combination x — vt and 
therefore dt = —vdx- The polar coordinate system that 
moves with the crack tip is related to the rest frame by 
r = y(x^^utp"+y2 and 9 = ta.n^'^[y/{x-vt)]. Thus, the 
traction-free boundary conditions on the crack faces are 
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(2) 



To proceed, we note that in the measurement region 
of [J] the maximal strain levels are 0.2 — 0.35 (see be- 
low) as the velocity of propagation varied from 0.20cs 
to 0.78cs, where Cs = \J pj p is the shear wave speed (/i is 
the shear modulus). These levels of strain motivate a per- 
turbative approach where quadratic elastic nonlinearities 
must be taken into account. Higher order nonlinearities 
are neglected below, though they most probably become 
relevant as the crack velocity increases. We write the 
displacement field as 



u{r,9) 



^^\r,9)^^u^'^\r,9)^0{^) , (3) 



where e quantifies the (dimensionless) magnitude of the 
strain. For a general U{F), s and <p can be expressed in 
terms of u of Eq. ([3]). Substituting these in Eqs. (H])-© 
one can perform a controlled expansion in orders of e. 



To make the derivation concrete, we need an explicit 
U{F) that corresponds to the experiments of [4|. The 
polymer gel used in these experiments is well-described 
by a plane stress incompressible Neo-Hookcan constitu- 
tive law [6|], defined by the energy functional [7,] 

C/(F) = |[F,,i^,,+det(F)-2-3] . (4) 

Using this exphcit U{F), we derive the first order 
problem in e 

^v^M^i) + SfiV(y ■ «(!)) = piiti' , (5) 

with the boundary conditions at 6 = ±'k 

r-^dgu^P+drU^^'^ = 0, Ar-^dBU^y^+2drU^^^ = . (6) 

This is a standard LEFM problem [3J. The near crack- 
tip (asymptotic) expansion of the steady state solution 
for Mode I symmetry is [3| 
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+ 0(r3/2).(7) 



Here Kj is the Mode I "stress intensity factor" and T is 
a constant known as the "T-stress" . Note that these pa- 
rameters cannot be determined by the asymptotic analy- 
sis as they depend on the global crack problem. fl{9; v) is 
a known universal function [3|, |8| . e in Eq. ([3]) can be now 
defined exphcitly as e = Kj/[4:fi^2Tri{v)], where i{v) is 
a velocity-dependent length-scale. i{v) defines the scale 
where only the order e and e^ problems are relevant. It is 
a dynamic length-scale that marks the onset of deviations 
from a linear elastic constitutive behavior. 

The solution of the order e equation, i.e. Eqs. ([7]), can 
be now used to derive the second order problem in e. The 
form of the second order problem for an incompressible 
material is 



^iV^u'-'^^ +3nV{V-u 



(2)N , M^g(^; v) _ ••(2) 



pu^ 



(8) 



The boundary conditions at 9 — zLtt become 



"^ (9) 
where contributions proportional to T were neglected. 
Here g{9;v) and k{v) are known functions, see [8\ and 
below. 



The problem posed by Eqs. ©-(121) has the structure of 
an effective LEFM problem with a body force ocr~^ and a 
crack face force ex r^^. Note that Eqs. ([S])-® are valid in 
the range ~£(u), where e^ is non-negligible with respect 
to e, but higher order contributions are negligible. Since 
one cannot extrapolate the equations to smaller length- 
scales, no real divergent behavior in the r — > limit is 
implied. We stress that the structure of this problem is 
universal. Only g{9; v) and k(v) depend on the second 
order elastic constants resulting from expanding a gen- 
eral U{F) to second order in e. For example, the oc r~^ 
effective body-force in Eq. ([5]) results from terms of the 
form d{du^^'du^^'), which are generic quadratic nonlin- 
earities. 



We now focus on solving Eq. ([5]) with the boundary 
conditions of Eqs. (O for the explicit g{9;v) and k{v) 
derived from Eq. ((4]) [sf. Our strategy is to look for a 
particular solution of the inhomogeneous Eq. ([5]) without 
satisfying the boundary conditions of Eqs. Q and then 
to add to it a solution of the corresponding homogeneous 
equation that makes the overall solution consistent with 
the boundary conditions. We find that the inhomoge- 
neous solution, T(6';w), is r-independent. The homoge- 
neous solution is obtained using a standard approach [3| 
by noting that the second boundary condition of Eqs. ([9|) 
requires that its first spatial derivative scales as r~^. The 
solution of the second order problem for Mode I symme- 
try is 



e'u^^\r,9;v) 
e^u'^\r,9-v) 
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Aad9d-B9s + Ty{9;v)^, t&n9d,s ^ ad.sta.n9, a^_, = 1 - wV^, , (10) 
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where A= [2asi? — 496)Tj,(7r; w) — K(-y)]/(2 — 4a^) (cf. Eq. The analytic form of T{9]v) depends mainly on 

([5])) and Cd is the dilatational wave speed. Cd = '2-Cs in an 
incompressible material under plane stress conditions. 



g{9\ v). The latter can be represented as 



9x{d]v)'. 



N(v) 
n=l 



an{v)cos{n9), gy{9;v) 



N(v) 
n=l 



bn(v) sin(ri0). 

(11) 

For u = we have -/V(O) = 3 and the representation is 
exact, while for higher velocities it provides analytic ap- 
proximations with whatever accuracy needed. For v ~ 
0.8cs only seven terms provide a representation that can 
be regarded exact for any practical purpose [8|. T{9;v) 
is then obtained in the form 



N(v) 



T^{9;v)~^Cn{v)cos{ne), Ty{e;v): 



n=l 



N{v) 



(u)sin(n0), 

(12) 

where the unknown coefhcients are determined by solving 
a linear set of equations [3|. A striking feature of Eqs. 
(fTO| is that they lead to strain contributions that vary 
as r~^, which are "more-singular" than the r~^'^ strains 
predicted by LEFM. 

We now show that the second order solution of Eqs. 
(fTO| entirely resolves the discrepancies raised by trying 
to interpret the experimental data of [J] in the frame- 
work of LEFM. The complete second order asymptotic 
solution, Eqs. ([3]), ([7]) and pO)) . contains three param- 
eters {Ki, T and B) that cannot be determined from 
the asymptotic solution and therefore must be extracted 
from the experimental data. 

These parameters were chosen such that Eqs. ([3]), ([7]) 
and (llOp properly describe the measured Ux{t, 0). Exam- 
ples for v/cs = 0.20, 0.53 and 0.78 are provided in Fig. 
[1] (top). With Kj, T and B at hand, we can now test 
the theory's predictions for eyy{r,Q) with no adjustable 
free parameters. The corresponding results are compared 
with both the measured data [J] and LEFM predictions in 
Fig. [1] (bottom) . In general, the agreement with the ex- 
perimental data is excellent. These results demonstrate 
the importance of the predicted r~^ strain terms near the 
crack tip. i is estimated as the scale where the largest 
strain component reaches values of 0.10 — 0.15. For the 
data presented in Fig. [T^, Syy > Sxx, where Sxx — dxUx 
is obtained by differentiating Ux- Thus, (, can be read 
off of the bottom panel to be ~ 0.5 — 1mm. Similar es- 
timates can be obtained for every w, though not always 
does eyy>exx, e.g. Fig. [It. 

For V = 0.53cs (Fig. [T)d) the theory still agrees well 
with the measurements, although some deviations near 
the tip are observed. These deviations signal that higher 
order corrections may be needed, though second order 
nonlinearities still seem to provide the dominant correc- 
tion to LEFM. For higher velocities, it is not clear, a- 
priori, that second order nonlinearities are sufficient to 
describe the data. In fact, the strain component exxii^, 0) 
for 11 = 0.7803 reaches a value of '--^ 0.35 in Fig. [T]:, suggest- 
ing that higher order nonlinearities may be important. 




FIG. 1: (Color online) Top: Measured Ux{r,0) (circles) fitted 
to the X component of Eq. @ (solid line) for (a) v — 0.20cs 
with Ki = 1070PaV7n, T = -3150Pa and B = 18. (b) v = 
0.53cs with A"/ = 1250Paym, r=-6200Pa and 5 = 7.3 and 
(c) w = 0.78c, with Ki= 980Pa^, r = -6900Pa and B = 26. 
Bottom: corresponding measurements of eyy{r, 0) = dyUy{r, 0) 
(circles) compared to the theoretical nonlinear solution (cf. 
Eq. (|3]l) with no adjustable parameters (solid lines); Ki, T 
and B are taken from the fit of Ux{r, 0). (dashed lines) LEFM 
predictions (analysis as in [J]) were added for comparison. 



Nevertheless, the second order theory avoids a funda- 
mental failure of LEFM; at high velocities {v>0.73cs for 
an incompressible material) LEFM predicts (dashed line 
in Fig. [TJ;) that the contribution proportional to Kj in 
eyy{r,0) (derived from Eqs. (O) becomes negative. This 
implies that eyy{r,0) decreases as the crack tip is ap- 
proached and becomes compressive. This is surprising, 
as material points straddling y — must be separated 
from one another to precipitate fracture. Thus, the sec- 
ond order nonlinear solution (solid line), though applied 
beyond its range of validity, already induces a qualitative 
change in the character of the strain. This is a striking 
manifestation of the breakdown of LEFM, demonstrating 
that elastic nonlinearities are generally unavoidable, es- 
pecially as high crack velocities are reached. The results 
of Figs. [T^-c both provide compelling evidence in favor of 
the developed theory and highlight inherent limitations 
of LEFM. We note that £{v) increases with increasing v, 
reaching values in the mm-scale at very high v. 

Our results indicate that the widely accepted assump- 
tion of "K-dominance" of LEFM, i.e. that there is always 
a region where the r"^ " strain term dominates all other 
contributions, is violated here. The results presented in 
Fig. [T] explicitly demonstrate that quadratic nonlineari- 
ties become important in the same region where a non- 
negligible T-stress exists. As elastic nonlinearities inter- 
vene before the r~^'^ term dominates the strain fields, 
the contributions of both of these terms must be taken 
into account as one approaches the crack tip. Since val- 
ues of the T-stress and of B are system specific, this 
observation is valid for the specific experimental system 
under study. They do indicate that the assumption of 
"K-dominance" is not always valid. 

An additional puzzle raised in [J] was that although 
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FIG. 2: (Color online) Measured crack tip profiles ((^j;(r, ±7r) 
vs. (jf>i(r, tt)) (circles). Shown are the parabolic LEFM best fit 
(dashed line) and the profiles predicted by the second order 
nonlinear corrections (solid line), (a) u = 0.2ca and (b) v = 
O.SSCs. T and B are as in Fig. [T] In contrast to the ~ 20% 
discrepancy in values of Ki obtained in [j), the respective 
values Ki — 1170Pay'm and Ki — ISOOPay^ correspond to 
within 9% and 4%, respectively, of Ki obtained from Ux(r,0) 
using the nonlinear theory, cf. Fig. [l] 



the form of both Ux {r, 0) and the Crack Tip Opening 
Displacement (CTOD) agreed with LEFM, the respec- 
tive derived values of Kj differed by about 20%, cf. Fig. 
3a in 4]. This puzzle is resolved by the theory as follows. 
The form of the CTOD is given by (py (r, ±7r) as a function 
of the distance, (j)x{''', t''), from the crack tip in the moving 
(laboratory) frame. Substituting 9 = tt into Eqs. ^, ([7]) 
and (fTU|) . the nonlinear theory predicts that the CTOD 
remains parabolic, where the log{r) term in (j)x{r,Tr) is 
negligible compared to r. This occurs at the same scale 
£{v) at which nonlinear corrections are essential to de- 
scribe the strain at 9 = 0, cf. Fig. 1. Quantitatively, the 
parabolic CTOD can be described with Kj values that 
differ from those describing Ux{r,0) by only a few per- 
cent with the same values of T and B (cf. Fig. [T]). This 
small Kj variation is possibly related to sub-leading non- 
linear corrections associated with the T-stress and will be 
addressed elsewhere. 

Let us now consider the CTOD in the near vicinity 
of the crack tip, i.e. when r is further reduced. Eqs. 
(fTO| predict the existence of log-terms in (px {r, 0) ■ These 
terms, which are negligible at = 7r on a scale i?(w), must 
become noticeable at smaller scales. Although this region 
is formally beyond the range of validity of the expansion 
of Eq. ([3]) , we would still expect the existence of a CTOD 
contribution proportional to logr to be observable. We 
test this prediction in Fig. [2] by comparing the mea- 
sured small-scale CTOD to both the parabolic LEFM 
form and the second order nonlinear solution with no ad- 
justable parameters. We find that these log-terms, whose 
coefficients were determined at a scale (.{v), capture the 
initial deviation from the parabolic CTOD at 6' = ±7r to 
a surprising degree of accuracy. This result lends further 
independent support to the validity of Eqs. (|10[) . 

In summary, we have shown that the second order so- 



lution presented in Eqs. pO|) resolves in a self-consistent 
way all of the puzzles that were highlighted in [J] . This 
solution is universal in the sense that its generic prop- 
erties are independent of geometry, loading conditions 
and material parameters. We would entirely expect that 
any material subjected to the enormous deformations 
that surround the tip of a crack must experience at least 
quadratic elastic nonlinearities, prior to the onset of the 
irreversible deformation that leads to failure. Our results 
show that these deformations, which are the vehicle for 
transmitting breaking stresses to crack tips, must be sig- 
nificantly different from the LEFM description, especially 
at high V. 

One may ask why we should not consider still higher 
order elastic nonlinearities. We surmise that quadratic 
elastic nonlinearities may be special, as they mark the 
emergence of a dynamic length-scale l{v) that character- 
izes a region where material properties - like local wave 
speeds, local response times and anisotropy - become de- 
formation dependent. This line of thought seems con- 
sistent with the observations of Refs. [9J. As support- 
ing evidence for this view, we note that the geometry- 
independent wave-length of crack path oscillations dis- 
cussed in [2, llO| seems to correlate with the mm-scale 
(.{v) at high V. Therefore, our results may have implica- 
tions for understanding crack tip instabilities. 
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